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We study bar formation in galactic disks as a consequence of the collision of two spiral galaxies 
under the influence of a potential which is obtained from the Newtonian limit of a scalar-tensor 
theory of gravity. We found that dynamical effects depend on parameters (a, A) of the theory. 
In particular, we observe that the bar is shorter for weaker tidal perturbations, which in turn 
corresponds to smaller values of A used in our numerical experiments. 

PACS numbers: 0.50.+h: 98.62.Hr; 98.65.Fz 
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I. INTRODUCTION 

Observations of spiral galaxies indicate that the pres- 
ence of a central structure called a bar is a common 
featured!]. The instability of isolated stellar disks in 
galactic models leads to bar formation and is character- 
ized by Toomre's stability parameter QQ. The models 
with Q < 1 are subject to bar formation. However, we 
arc interested in dynamical effects of non isolated sys- 
tems which are found in clusters of galaxies. In this 
regard, it has been suggested that the observed bar in 
many spirals is the result of the gravitational interac- 
tion between two or more nearby galaxies. For instance, 
Nogushi[|[ has found that, during the collision of two 
galaxies and between the first and the second closest ap- 
proaches, the disk takes on transient bar shape. The 
gravitational interaction between the two galaxies gives 
rise to perturbations in the orbits of the stars that results 
in the formation of the bar. 

Bar formation in the simulations of stellar disks de- 
pends upon various simultaneous effects. In the case 
of collisions, these factors areQ: rotation curve shape, 
disk-halo mass ratio, perturbation force and geometry. 
However, simulations can suffer from numerical effects 
such as low spatial and temporal resolution, too few par- 
ticles representing the system and an approximate force 
model. These effects can be drastic: for example, obser- 
vations show that bars have typically a length scale close 
to the exponential length of the disk[5(, while the bar's 
semi-major axis obtained from numerical models is two 
to four times longer [E0] ■ 

Recent observational data measured in the Cosmic Mi- 
crowave Background at various angular scales, in Super- 
novae la, in the 2dF galactic survey, and baryon acoustic 
oscilations, suggest [8, 9] that the Universe is composed of 
about 4.6% baryons in the form of gas and stars, 23.2% 
dark matter (DM) and 72.2% dark energy, which is a 
kind of cosmological constant and is responsible for the 
today accelerated expansion of the Universe. In this way, 
galaxies are expected to possess these dark components 
and, in accordance with rotation curves of stars and gas 



around the centers of spirals, this might be in the form 
of halos, and must contribute to at least 3 to 10 times 
the mass of the visible matter of spirals. 



Regarding to the nature of DM, we know that DM 
has to be non-baryonic. This is because nucleosynthesis 
abundances of light elements are only consistent with the 
above-mentioned baryonic fraction, and this is not suffi- 
cient at all to account for rotational velocities of spirals. 
This fact opens up new possibilities for explaining the 
nature of DM. In this sense, in papers fiol . fill . H2I [l3| a 
model was proposed in which a scalar field (SF) couples 
non-minimally to gravity to produce locally a modified 
Newtonian theory of gravity. It turns out that the dy- 
namics is now determined by the Poisson equation cou- 
pled to a Klein-Gordon equation for the assumed scalar 
field in the galaxy. Thus, the boson mass of the scalar 
field modifies the Newtonian law of attraction, and the 
dynamics of DM is different from its Newtonian counter- 
part. In this scalar-tensor theory potential-density pairs 
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for various halo density profiles were computed 



In the present work we use the above-mentioned results 
to study the collision process of two spirals, each of which 
possess a disk, bulge and dark halo, in order to estimate 
the effects of the modified gravity theory on the bar's 
length and the orbital decay of galaxies. A numerical 
treecode was developed by one of us (MARM) in order 
to compute the evolution of an A-body system interact- 
ing with the standard gra vity plus the interaction of the 
scalar field [IE EI EE EE Gj • We first study dynamical 
effects on isolated galaxy models for three different inter- 
action scales (A). We found no significant changes in the 
morphology of models, but we did in the total potential 
energy. Then, we analyze the formation of a bar during a 
parabolic collision of two identical galaxies and compare 
the results obtained in the SF model with three scales of 
A and with a pure Newtonian interaction. 
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II. SCALAR-TENSOR THEORY AND ITS 
NEWTONIAN LIMIT 

A typical scalar-tensor theory is given by the following 
Lagrangian 



C = 



16tt 



-<t>R- 



(d<f>) 2 - V{ct>) 



+^m(9hu) , (1) 



from which we get the gravity and SF equations. Here 
<7 M „ is the metric, Cm{.9iiv) is the Lagrangian matter and 
u>(4>) and V((j>) are arbitrary functions of the SF. Ac- 
cording to the Newtonian approximation, gravity and 
SF are weak, and the velocities of the stars are non- 
relativistic. Then, we expect to have small deviations of 
the SF around the background defined here as (</>) = Gq 1 . 
If we define the perturbation = <fi — (0) , then the New- 
tonian approximation gives the equations [l2l [l8| 



\7 2 ip = AttGqp , \J 2 <j) - m 2 <f) = -8naG p 



(2) 



where tp = 4(^00 + </>)• Here we define A = K/mc, the 
Compton wavelength of the effective mass m of some el- 
ementary particle (boson) given through ui{4>) an d the 
potential V(cf>), and a = 1/(3 + 2lu(4>)) is the amplitude 
of the perturbed SF, <j). The above formalism is valid 
for any potential that can be expanded in Taylor series 
around (<j>). In what follows we will use A instead of 
?7i _1 . This mass can have a range of values depending on 
particular particle physics models. 

General solutions of Eqs. ([2]) can be found in terms of 
the corresponding Green's functions and the new Newto- 
nian potential is 



N 



1>-\4> = -G Q 



dr., 



- aGo I dr s ^P : +B.C. 



Solutions of these equations for point masses are 



A' 



where 



= —u — au\ , 
G m s 



E 



Ux = ^^£!!^ e -|r-r s |/A 



r - r, 



(3) 



(4) 
(5) 



(0) 
(7) 



with m s being a source mass. The total gravitational 
force on a particle of mass rrii is 



F = -V$at = mi a. 



(8) 



The potential u is the Newtonian part and u\ is the SF 
modification which is of Yukawa type. 



III. INITIAL CONDITIONS 

We use the Monte-Carlo procedure to construct a 
galaxy model with a Newtonian potential. A fully self- 
consistent model in the context of the SF is in prepa- 
ration. The initial conditions of the galaxies are con- 
structed following the model described by Barnes [nj. In 
this model, both the bulge and halo are non-rotating, 
spherically symmetric and with an isotropic Gaussian 
distribution of velocities characterized by the velocity 
dispersions o\, and er^, respectively. The units are such 
that the local (r -C A) gravitational constant is G = 
Gq(1 + a) — 1, and the units of mass, longitude and time 
are M = 2.2 x 10 11 M , R = 40 kpc and T = 250 Myrs, 
respectively. The bulge density profile is[2fj| 



Pb(r) 



M b a b 



1 



27r r(r + ab) 3 



(9) 



and the halo density profile is a Dehnen's family member 
with 7 = [HI 



Ph(r) 



3M h 



Air (r + a h ) 4 
The disk density profile is exponential [2z 



Pd(r,z) 



M, 



- r / Qd sech 2 f — 



(10) 



(11) 



Here M b = 0.0625, M d = 0.1875, and M h = 1.0 are 
the total mass of the bulge, disk, and halo, respec- 
tively. The scale lengths of the bulge, halo, and disk 
are a& = 0.04168, ah — 0.1, and a<j = 1/12, respec- 
tively, and zq — 0.007 is the scale height of the disk. The 
mass distributions were truncated at a radius contain- 
ing 95% of the total mass, since they extend to infin- 
ity. The compound galaxy was sampled with N — 40960 
equal mass particles. The velocity distribution of the 
disk is given by the Schwarzschild distribution function 
with the velocity dispersions = 2a z oc e~ adT , and a z 
given by the equilibrium condition of an infinite gravitat- 
ing sheet; a<f, is calculated from the epicyclic approxima- 
tion. The Toomre's parameter of initial disk is Q « 1, so 
we have the disk which is marginally stable for axisym- 
metric perturbations, but not, however, against strong 
non-axisymmetric ones. 

Observations suggest that the majority of the interact- 
ing galaxies are located on nearly parabolic orbits. For all 
the collisions, disks were located in the plane of parabolic 
orbits, calculated from parameterized equations of the 
two-body problem, with a pericentric separation p = 0.4, 
and the time to pericenter t p = 3.0. The direction of 
rotation of one of the disks (disk 1) was in the same di- 
rection as the corresponding orbital-angular-momentum, 
i.e., direct motion. The other disk (disk 2) was in ret- 
rograde motion. The two colliding galaxies are initially 
identical. 
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IV. NUMERICAL METHOD 

For the time evolution we use a Barnes tree code 
type Hl| modified to include the Newtonian contribu- 
tion of the scalar fields as given by Eqs. (4)-(8) [ilj . 
The forces were computed with a tolerance parameter 
8 = 0.75, and including the monopole term only. For the 
gravitational potential we used the standard Plummer 
model 



$ oc -- 



1 



(12) 



Here, e is the softening parameter taken in our simula- 
tions to be e = 0.015. The equations of motion were in- 
tegrated using the second order leap-frog algorithm with 
a fixed time step At = 1/256. With these parameters we 
obtain a good energy conservation (< 0.2%) and, also, 
good angular momentum conservation (< 0.5%) for all 
runs presented here. 

To characterize quantitatively the bar amplitude, we 
consider the distortion parameter defined as 24 



where 



and 



V+ 



Ixx lyy 



N 



2h 



(13) 



(14) 



(15) 



fe=i 



The particles that are outside of the spatial region of the 
original disk can affect the parameter under study. For 
instance, if we calculate the distortion parameter using 
all the particles in the disk, we have in both cases sim- 
ilar evolution curves. Therefore, to avoid the noise, we 
exclude particles that are outside of the original radius 
of the disk. 



V. RESULTS 

We first study isolated galaxy models with different 
values of A. We consider four set of simulations followed 
up to time t = 8.0. The parameters and results of runs 
are presented in the table, where Eq and E are the ini- 
tial and mean total energies, and AE/Eq is the relative 
change of the total energy during the evolution with re- 
spect to its initial value. Though the scale of interaction 
A and magnitude of a are unknown, we choose their val- 
ues arbitrarily, such that A is equal to the cutoff radius 
of the disk, bulge and halo for models Al, A2 and A3, 
respectively, and a fixed amplitude of the SF, a = 1. The 
larger A makes weaker the contribution of SF for a fixed 
galactic size. For A = oo one obtains the Newtonian case, 



model A4. Previous studies of protogalactic interactions 
under influence of this SF[10| were made for scales less 
than those considered here. 



Run 


A 


Q 


\Eo\ 


\E\ 


i w i. % 


V 


Al 


0.4 


1.0 


0.7279 


0.7277 


0.066 


0.030 


A2 


1.0 


0.9 


0.9511 


0.9509 


0.048 


0.033 


A3 


6.0 


0.9 


1.1632 


1.1630 


0.039 


0.042 


A4 


oo 


1.0 


1.2234 


1.2232 


0.036 


0.045 



All models show a good energy conservation (see table). 
The presence of the SF decreases the total potential en- 
ergy due to shallower potential well at distances r > A. 
The initial models reaccommodate rapidly due to poten- 
tial modification, i.e., shifts to a new equilibrium state. 
At the end of evolution, the components of galaxy mod- 
els with the SF became slightly more extended. The 
velocity profiles of the components match the Newtonian 
ones up to r « A. For r > A there is a slow decay in 
velocities, since the effective gravitational constant de- 
creases. The distortion parameter shows a nearly equal 
noise level of surface density of the disks. The evolution 
of the Toomre's local stability parameter Q shows a slow 
decay from Q rj 1 to Q rs 0.9 with mean values presented 
in the table for each run. 

Then, we proceed to study the interaction of two equal 
galaxy models. During the orbital decay we analyze the 
bar's strength for different values of A. Because the equi- 
librium galaxy models were constructed with Newtonian 
potential, we relax them up to time t = 1.0 with modified 
SF potential in order to reach a new equilibrium state for 
a given A. Then we place relaxed galaxies on parabolic 
orbits and let them interact. 

The results are as follows. The first encounter occurs at 
time t w 3.0. This is the time of major transfer of orbital 
angular momentum, and then disk 1 forms a strong bar, 
while disk 2 is in retrograde motion and develops a very 
weak bar. Figure 1 (left) presents the separation between 
the centers of mass, p, of two galaxies as a function of 
time. All runs, except Al, after several close approaches 
merge and form a remnant. In the run with model Al, 
after the first approach the galaxies separate such a large 
distance that they will never encounter again. This is be- 
cause the weaker gravity is, the more lessable the galax- 
ies are to become bound: at distances larger than A the 
potential diminishes in comparison with the pure New- 
tonian. In the runs A2 and A3 the SF causes just a 
retardation of the subsequent interactions. The simula- 
tion with model A3 is practically identical to Newtonian 
one. 

In order to analyze the bar formation, we consider fol- 
lowing only the disk 1. In figure 1 (right) we plot the 
distortion parameters r\ as a function of time for runs 
with galaxy models A2 and A4 only. As a consequence 
of gravity modification, the galaxies do not approach each 
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other too closely as happens in the Newtonian run. Thus, 
weaker perturbations make the bar shorter. The run with 



VI. CONCLUSIONS 

From simulations of isolated galaxy models with differ- 
ent A, we can see that the addition of a non-minimally 
coupled SF slightly modifies the equilibrium of Newto- 
nian model, acting as a small perturbation, and it dimin- 
ishes the total potential energy for r > A, since the effec- 
tive gravitational constant decreases in this range. Our 
results show that the interaction of galaxies with the SF is 
weaker in comparison with the Newtonian case. We have 
found that the inclusion of the SF changes the dynami- 
cal properties such as the collision time, bar morphology, 
and in general the remnant properties. All these changes 
depend on the pair (a, A), which on the other hand, can 
be constrained from observations. For instance, the du- 
ration of interaction cannot be larger than the age of the 



galaxy model A3 is very similar to the Newtonian case. 
The bar's phases are displaced due to orbit modification. 

I 

Universe, implying constraints on values of Go, which 
depends on G and a. These constraints can be provided 
from statistical data on the fraction of observed inter- 
acting galaxies. A wide range of parameters should be 
investigated and higher resolution have to be used in sim- 
ulations in order to make predictions for particular inter- 
acting models. Further investigations with more particles 
and self-consistent initial models are under way. 
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